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ABSTRACT 

A  few  years  ago,  Lee  and  Carder1  demonstrated  that  for  the  quantitative  derivation  of  major  properties  in  an 
aqua-environment  (information  of  phytoplankton  biomass,  colored  dissolved  organic  matter,  and  bottom 
status,  for  instance)  from  remote  sensing  of  its  color,  a  sensor  with  roughly  ~17  spectral  bands  in  the  400  - 
800  nm  range  can  provide  acceptable  results  compared  to  a  sensor  with  81  consecutive  bands  (in  a  5-nm 
step).  In  that  study,  however,  it  did  not  show  where  the  17  bands  should  be  placed.  Here,  from  nearly  300 
hyperspectral  measurements  of  water  reflectance  taken  in  both  coastal  and  oceanic  waters  that  covering  both 
optically  deep  and  optically  shallow  waters,  first  and  second  derivatives  were  calculated  after  interpolating 
the  measurements  into  1-nm  resolution.  From  these  hyperspectral  derivatives,  the  occurrence  of  zero  value  at 
each  wavelength  was  accounted  for,  and  a  spectrum  of  the  total  occurrences  was  obtained,  and  further  the 
wavelengths  that  captured  most  number  of  zeros  were  identified.  Because  these  spectral  locations  indicate 
extremum  (a  local  maximum  or  minimum)  of  the  reflectance  spectrum  or  inflections  of  the  spectral 
curvature,  placing  the  bands  of  a  sensor  at  these  wavelengths  maximize  the  possibility  of  capturing  (and  then 
accurately  restoring)  the  detailed  curve  of  a  reflectance  spectrum,  and  thus  maximize  the  potential  of 
detecting  the  changes  of  water  and/or  bottom  properties  of  various  aqua  environments  with  a  multi-band 
sensor. 


Keywords:  Ocean-color  remote  sensing,  spectral  bands 
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1.  INTRODUCTION 

Since  the  successful  demonstration  of  the  Coastal  Zone  Color  Scanner  (CZCS)  in  measuring  the  spatial  and 
temporal  variation  of  ocean  color2, 3,  the  importance  of  observing  ocean  and  coastal  waters  with  sensors  in 
the  visible  domain  is  getting  more  and  more  attention  from  various  countries.  Especially,  after  the  launch  of 
the  SeaWiFS  and  MODIS  satellites  about  ten  years  ago  and  the  successful  and  exceptional  results  of  such 
data  since  then,  more  and  more  satellite  sensors  have  been  (or  planned  to  be)  launched  for  the  observation  of 
global  ocean  via  observing  the  water  color.  These  include  the  COCTS  of  China,  the  MERIS  of  ESA,  the 
OCM  of  India,  and  the  GOCI  of  Korea4,  etc. 

Ideally,  to  maximize  the  potential  of  observing  properties  of  aqua-environment  and  their  temporal  and 
spatial  changes,  it  is  best  to  have  a  sensor  capable  of  colleting  signals  with  a  continuous  spectrum.  After 
consideration  of  data  flow,  processing,  and  storage  as  well  as  signal-to-noise  ratio,  satellite  sensors 
commonly  have  a  few  bands  in  the  visible  domain4.  For  instance,  the  CZCS  had  five  bands,  while  SeaWiFS 
and  MODIS  have  six  bands  (with  slightly  different  configurations).  Selection  of  those  bands  (number  of 
bands  and  their  center  wavelengths  as  well  as  band  widths)  was  based  on  the  optical  characteristics  (both 
absorption  and  scattering)  of  phytoplankton  and  colored  dissolved  organic  matter4.  Thus  these  configurations 
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are  in  general  good  for  monitoring  chlorophyll  concentrations  (index  of  phytoplankton)  and  CDOM  in  the 
open  oceans,  but  not  optimal  for  remote-sensing  of  complex  coastal  waters.  Lee  and  Carder  demonstrated 
that  for  adequate  remote  sensing  of  major  properties  (water  column  and  bottom)  of  both  open  ocean  and 
complex  coastal  aqua-environments,  a  sensor  requires  ~  17  spectral  bands  in  the  400  -  800  nm  range  to 
obtain  similar  results  as  a  sensor  with  81  consecutive  bands  in  a  5-nm  spacing.  It  is  not  yet  clear,  however, 
where  the  1 7  bands  should  be  placed. 

In  this  study,  with  an  aim  for  future  advanced  sensor  for  observing  properties  of  most  aqua-environments 
from  water  color,  we  derived  the  spectral  locations  of  the  optimal  bands  through  first-  and  second  derivatives 
of  the  remote  sensing  reflectance.  These  bands  thus  can  capture  most  of  the  spectral  variations  of  remote¬ 
sensing  reflectance  caused  by  varying  water  column  and/or  bottom  properties. 

2.  DATA  AND  METHOD 

Since  1993,  in  a  life  span  of  more  than  12  years,  an  extensive  data  set  of  hyperspectral  (-350  -  900  nm  with 
a  resolution  -2  nm)  remote-sensing  reflectance  for  various  types  of  water  was  collected,  with  most  of  them 
measured  by  the  group  lead  by  Prof.  Kendall  Carder  at  the  College  of  Marine  Science  of  the  University  of 
South  Florida.  Table  1  summarizes  the  cruises  and  general  information  of  the  various  environments.  The 
measured  waters  included  open  ocean  blue  waters  to  turbid  coastal  yellowish  waters  and  optically  shallow 
waters.  Wide  ranges  of  concentrations  of  chlorophyll  and  dissolved  organic  matter  and  suspended  sediments, 
along  with  different  types  of  phytoplankton  taxonomy  were  encountered. 


Table  1:  Data  used  in  this  study  (a  total  of  298  stations) 


AREA 

DATE 

#  OF  STATIONS 

rCHLl  RANGE 

Gulf  of  Mexico 

Apr  1993 

24 

0.07-49.0 

Florida  Keys 

July  1994 

5 

0.06  -  0.5 

Arabian  Sea 

Dec  1994 

20 

0.3 -0.9 

Chesapeake  Bay 

Sept  1996 

36 

1.7-20.7 

Hawaii 

Feb  1997 

6 

0.1 -0.3 

Florida  Keys 

May  1997 

8 

N/A 

East  China  Sea 

July  1998 

37 

0.5 -2.8 

California  Coast 

Oct  1999 

37 

0.2  -  9.4 

North  Atlantic 

July  2001 

17 

4.5-13.9 

Monterey  Bay 

Apr  2003 

56 

0.1 -9.7 

|  Ft.  Lauderdale 

July  2005 

52 

N/A 

Remote-sensing  reflectance  (Rrs(k))  is  defined  as  the  ratio  of  water-leaving  radiance  to  downwelling 
irradiance.  For  each  station  of  the  many  cruises,  using  a  handheld  spectroradiometer  (Spectron  Model  SE- 
590  of  Spectron  Engineering  before  1997;  SPECTRIX  for  1997  and  later),  a  series  (—  3  to  5  scans)  of 
upwelling  radiance  above  the  surface  (Lu(0+))  and  the  downwelling  sky  radiance  (Lsky)  were  directly 
measured  by  the  instruments.  The  measured  Lu  is  a  sum  of  photons  emerged  from  subsurface  scattering  (the 
desired  signal)  plus  surface  reflected  sky  and  solar  light  (noise).  Downwelling  irradiance  was  derived  by 
measuring  the  radiance  (JLr)  reflected  from  a  standard  diffuse  reflector  (Spectralon®).  For  each  of  the 
collected  scans,  hyperspectral  total-remote-sensing  reflectance  (Tra(X))  and  sky  reflectance  (5„(X.))  were 
derived  through 
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(1) 


L  Rn  , 
Trs  =  — — -  and 

Lr  * 

with  Rr  the  reflectance  of  the  diffuse  reflector  (~  10%). 


^sky  Rr 

Lr  *  ’ 


Based  on  these  measured  Tn  and  SK  curves,  averages  of  7^  and  Sn  for  each  station  were  obtained, 
respectively,  in  order  to  reduce  the  random  variations  associated  with  measurements  due  to  reflection  from 
different  wave  facets,  etc.  Further,  because  7^  is  the  sum  of  water-leaving  radiance  and  radiance  reflected  by 
the  water  surface,  Rt s(X)  is  calculated  as5, 6 


R„W  =  T„(A)-FS„V l)-A,  (2) 

where  F  is  surface  Fresnel  reflectance  (around  0.023  for  the  viewing  geometry),  and  A  accounts  for  the 
residual  surface  contribution  (glint,  etc.).  Figure  1  shows  some  examples  of  spectral  Trs  and  Srs,  while  Figure 
2  shows  examples  of  Rn(X).  Clearly,  as  already  well  known,  Srs  spectra  are  monotonic  decreasing  functions 
of  X,  where  as  Trs  spectra,  because  of  the  contribution  of  RK,  show  distinct  variations  from  water  to  water. 
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Fig.  1.  Examples  of  Tn  and  Sn  spectra  used  in  this  study. 


Rn  is  a  function  of  the  total  absorption  and  backscattering  coefficients7'10.  Spectra  of  backscattering 
coefficients,  although  normally  stronger  in  the  shorter  wavelengths,  generally  do  not  have  strong  spectral 
signatures3,  11 .  The  spectrum  of  total  absorption  coefficient,  however,  is  an  additive  result  from  the 


contributions  of  CDOM,  various  pigments  contained  in 
phytoplankton,  plus  water  molecules3,  l2.  These 
contributions  are  spectrally  selective,  thus  resulted 
significant  variations  in  the  spectral  signatures  of  total 
absorption  coefficient13, 14.  The  extremum  values  (local 
maximum  or  minimum)  of  Rn(k)  and  the  inflections  of 
7?rs(X)  indicate  the  different  combinations  of  those 
optically  active  constituents  and/or  bottom  properties. 

If  a  sensor  has  only  a  few  spectral  bands,  to  best 
capture  the  spectral  signatures  of  Z?rs(X),  then  these 
bands  need  to  be  positioned  at  wavelengths  that  can 
match  the  locations  of  those  extremums  and  inflections 
of  /?rs(X).  The  spectral  locations  of  the  extremums  and 
inflections  of  /?re(X)  are  those  wavelengths  with  zero 


wavelength  [nm] 


Fig.  2.  Examples  of  used  in  this  study. 
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first-  and  second-order  derivatives  of  RIS(X),  respectively. 

To  find  out  the  optimal  spectral  locations  for  all  collected  samples,  we  carried  out  the  following 
calculations: 


a(A)  = 


dKM) 

d; i 
da(A) 
~dX 


(3) 


where  o(X)  and  q(X)  are  the  first-  and  second-order  derivatives  of  Rn(X),  respectively.  Before  calculating  the 
derivatives,  R^X)  spectrum  was  interpolated  into  1-nm  resolution  and  was  smoothed  with  a  5-nm  running 
average  (to  remove  some  of  the  random  noises).  Figure  3  presents  examples  of  a(X)  and  q(X.)  of  the  data  set. 


Fig.  3.  Examples  of  c(A.)x  10000  (Left)  and  q(X.)xl000  (Right). 


3.  RESULTS  AND  DISCUSSION 


From  all  of  the  calculated  o(X)  and  q(X)  spectra,  we  obtained  the  spectral  distribution  (in  1-nm  resolution) 
where  o(X)  and  q(X)  equal  zero,  respectively.  Mathematically,  these  distributions  are, 


° 

f°(»  =  ^Tr - 


f(W  = 


E<r,W  =  o 


(4) 


Here  cr,(X)  and  q,{X)  are  the  z'th  o(X)  and  q(X)  spectrum,  respectively.  Figure  4  shows  fo(k)  and  (after  a  5- 

nm  running  average  to  highlight  the  major  bands)  of  this  study. 

For  an  Rrs(X)  spectrum,  o(A.)  =  0  indicates  that  there  is  an  extremum  Rrs  value  at  X.  For  various  waters, 
because  of  the  different  combinations  of  water  constituents  and/or  bottom  characteristics,  the  locations  of  the 
extremums  shifts  (see  Fig.2).  The  spectral  distribution  of  fJX)  then  summaries,  for  all  data  in  this  study,  the 
wavelength  behavior  in  capturing  the  extremums.  The  greater  the/^X)  value,  the  more  the  appearance  of  an 


Proc.  of  SPIE  Vol.  6680  66800D-4 


extremum  at  that  X.  Apparently,  from  this  extensive  data  set,  the  bands  that  captured  more  extremums  of 
are  concentrated  to  seven  bands  (in  the  360  -  700  nm  range)  instead  of  scattered  to  all  bands  (see  Fig.4, 
Left).  The  seven  bands  are  centered  at  (round  to  nearest  0  or  5):  395,  440,  490,  515,  565,  665,  and  685  nm. 
Such  a  result  is  not  surprising.  For  instance,  440  nm  is  at  the  peak  absorption  of  phytoplankton  pigments 
(chlorophyll-a),  which  normally  causes  a  local  minimum  of  Rrs(X);  on  the  other  hand,  685  nm  is  around  the 
peak  of  chlorophyll  fluorescence,  which  normally  contributes  a  local  maximum  of  Rrs(X)  (see  Fig.2).  This 
result  from  the  analysis  of  the  first  derivative  is  consistent  with  the  results  of  the  principal  component 
analysis15  where  five  bands  in  the  400  -  600  nm  range  were  suggested  for  the  detection  of  phytoplankton. 
Therefore,  for  a  multi-band  sensor,  it  requires  that  it  has  at  least  these  bands  to  best  capture  the  extremums  of 
the  various  /?„( X)  spectra. 

q(X)  =  0,  on  the  other  hand,  indicates  an  inflection  of  an  R^X)  curvature  at  X,  normally  resulted  from 
different  combinations  of  pigments  and  water  molecules,  and/or  bottom  effects.  Table  2  summarizes  the 
wavelengths  that  capture  more  inflections  of  R^A.)  curves  (see  Fig.  4,  Right).  Different  from  the  distribution 
of  fa(X),  the  greater  values  of  f^X)  are  much  less  concentrated  spectrally.  Such  a  result,  however,  does 
indicate  that  our  data  set  covers  a  wide  range  of  water  environments,  which  encompass  various  types  of 
phytoplankton  classes  (resulted  minor  changes  in  R,S(X)  curvatures  because  of  pigment  absorption).  For  a 
multi-band  sensor,  it  is  necessary  to  have  bands  cover  these  wavelengths  if  pigment  compositions  (related  to 
phytoplankton  classes  and/or  functional  groups)  and/or  bottom  properties  are  also  desired  from  remote 
sensing  of  water  color.  Interestingly,  these  bands  match  very  well  with  the  center  wavelengths  of  the  major 
pigments  (Table  2)  of  presented  in  Hoepffner  and  Sathyendranath16. 


400  450  500  550  600  650  700 


400  450  500  550  600  650  700 


40 


Wavelength  [nm] 


Wavelength  [nm] 


Fig.  4.  Distributions  of /CT(X)  (Left)  and^(X)  (Right). 


Consolidating  the  major  bands  that  have  greater  values  of  fa(A.)  and  f?(X)  (and  considering  a  10-20  nm 
bandwidth  for  satellite  sensors),  possible  bands  for  future  sensors  are  suggested  (Table  2),  although 
absorption  effects  by  gases  in  the  atmosphere  are  not  considered  yet.  Interestingly,  the  analysis  from  the 
first-  and  second-order  derivative  also  indicates  a  number  of  16  bands  for  the  observation  of  ocean  color.  We 
added  a  band  at  710  nm,  consistent  with  the  configuration  of  MERIS,  although  not  a  big  number  of  fo(710) 
was  obtained  (because  of  low  number  of  observations  of  such  extreme  cases).  This  band,  however,  is 
important  for  turbid  coastal  water  and  water  with  extreme  phytoplankton  blooms  (red  tide)17' 18. 
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Table  2  also  presents,  as  an  example,  the  spectral  bands  of  MODIS  and  MERIS.  Clearly,  the  wavelengths 
of  those  bands  are  in  general  consistent  with  findings  presented  here.  Comparing  MODIS  bands  with  that 
from  this  study  and  that  of  MERIS,  the  significant  shortage  of  MODIS  (and  SeaWiFS)  is  a  band  between 
551  and  667  nm.  Such  a  band  is  of  great  value  for  remote  sensing  of  suspended  sediments19  or  remote 
sensing  of  optically  shallow  waters. 


Wavelength 

From  fa{X) 

From/;(X) 

Proposed 

Bands 

H&S 

Bands* 

MODIS 

MERIS 

385 

X 

1 

384 

395 

X 

2 

400 

X 

413 

412 

412 

425 

X 

3 

440 

X 

4 

435 

443 

443 

445 

X 

460 

X 

5 

464 

475 

X 

6 

490 

X 

7 

490 

488 

490 

510 

'  X 

8 

510 

515 

X 

520 

X 

532 

531 

545 

X 

9 

555 

X 

551 

565 

X 

10 

560 

580 

X 

11 

583 

615 

X 

12 

623 

620 

635 

X 

13 

645 

X 

644 

655 

X 

655 

665 

X 

14 

667 

665 

670 

X 

676 

685 

X 

15 

678 

681 

710 

16 

709 

750 

For  atmosphere  correction 

748 

754 

870 

For  atmosphere  correction 

870 

870 

4.  SUMMARY 

In  this  study,  based  on  an  extensive  measurements  of  remote-sensing  reflectance  of  various  aqua 
environments,  the  primary  bands  (in  the  -400  -  700  nm  range)  that  optimally  capture  the  spectral  signatures 
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of  spectral  Rrs  are  determined  via  first-  and  second-order  derivatives.  These  bands  are  in  general  cover  the 
operational  bands  of  MODIS  and  MERIS,  but  provide  suggestions  of  extra  bands  in  order  to  provide  more 
and  improved  results  for  remote  sensing  of  oceanic  and  coastal  waters. 
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